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Abstract 

The method used and results obtained are described for calculation of 
the cross section for the pion diffractive dissociation into two jets. The main 
new qualitative result is that the distribution of longitudinal momenta for 
jets is not simply proportional to the profile of the pion wave function ^(x), 
but depends on it in a complicated way. In particular, it is shown that 
under conditions of the E791 experiment, the momentum distribution of 
jets is similar in its shape for the asymptotic and CZ wave functions. It 
is concluded therefore that, unfortunately, the process considered is really 
weakly sensitive to the profile of the pion wave function, and the accuracy of 
data is insufficient to distinguish clearly between different models of <j)n(x). 

Comparison with the results of other papers on this subject is given. 
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1. The E791 experiment at Fermilab [1] has recently measured 
the cross section of the hard diffractive dissociation of the pion into two jets. 
In particular, the distribution of the total pion longitudinal momentum into 
fractions y\ and y 2 , (yi + y 2 ) = 1, between jets has been measured. The 
main purpose was to obtain in this way the information about the leading 
twist pion wave function ^(xi, x 2 ), which describes the distribition of quarks 
inside the pion in the longitudinal momentum fractions X\ and x 2 = 1 — X\. 

The hope was based on theoretical calculations of this cross section, see 
[2]- [5]. It has been obtained in these papers that the cross section is simply 
proportional to the pion wave function squared: da/dyi ~ | 1 2 - In such 
a case, it would be sufficient to measure only the gross features of da/dy 
to reveal the main characteristic properties of 4> n {%), and to discriminate 
between various available models of 4) n (x). 

Our main qualitative result [6]- [7] is that this is not the case. The real 
situation is much more complicated, with da/dy depending on cj) n (x) in a 
highly nontrivial way. We describe below in a short form our approach and 
the results obtained for this cross section. 

2. The kinematics of the process is shown in fig.l. We take 
the nucleon as a target, and the initial and final nucleons are substituted 
by two soft gluons with momenta q\ = (u + £)P and q 2 = {u — £)P, P = 
(P + P')/2, A = (qi — q 2 ). [] The lower blob in fig.l represents the generalized 
gluon distribution of the nucleon, G^(u) = G^(—u) [8]- [10]. 

The final quarks are on shell, carry the fractions y\ and y 2 of the ini- 
tial pion momentum, and their transverse momenta are: (kj_ + (qj_/2)) and 
(— + (q_i_/2)), q± <C k±, where q± is the small final transverse momentum 
of the target, while k± is large. The invariant mass of these two quark jets 
is: M 2 = k 2 ± /yiy 2 , and £ = M 2 /v, v = (s-u). 

The upper blob M in fig.l represents the hard kernel of the amplitude. 
According to the well developed approach to description of hard exclusive 
processes in QCD [11]- [14] (see [15] for a review), all hard gluon and quark 
lines in all diagrams have to be written down explicitly and substituted by 
their perturbative propagators. In other words, the hard momentum flow 
have to be made completely explicit and these hard lines of diagrams con- 
stitute the hard kernel M. They should not be hidden (if it is possible at 
all) as (a derivatives of) "the tails" of the unintegrated pion wave function 

1 The small skewedness, £ -C 1, is always implied. It is typically: £ ~ 10~ 2 , in the 
Fermilab experiment. 



2 



^tt(%, l±), or of the "unintegrated gluon distribution". This is, first of all, 
what differs our approach from previous calculations of this process in [4]- [5]. 

For calculation of M in the leading twist approximation and in the low- 
est order in a s , the massless pion can be substituted by two massless on 
shell quarks with the collinear momenta x\p n and X2Pn and with zero trans- 
verse momenta, as account of primordial virtualities and transverse momenta 
results only in higher twist corrections to M. The leading twist pion wave 
function (j) n (x,/j, ) describes the distribution of these quarks in momentum 
fractions X\ and x 2 - 

So, the hard kernel M is proportional to the scattering amplitude of two 
initial collinear and on shell quarks of the pion on the on shell gluon: 

M ~ {d(xxPn) + u(x 2 p 7T ) + g(qi) -> d(pi) + u(p 2 ) + g(q 2 )} • 

In lowest order in a s (k±) M consists of 31 connected Born diagrams, each 
one is ~ 0(a^(k±)). 

The general structure of the amplitude is, therefore (symbolically): 

T ~ (P'lA 1 - ■ A^P) ® ftJIftlij) ® (0|u ■ d\ir-), 

where the first matrix element introduces the skewed gluon distribution of 
the nucleon G^(u), ipi and ip 2 are the free spinors of final quarks, "M" is the 
hard kernel, i.e. the product of all vertices and hard propagators, the last 
matrix element introduces the pion wave function 4> n (x), and (g> means the 
appropriate convolution. 

3. Proceeding in the above described way and summing up con- 

2 As usual, on account of leading logs from loops the soft pion wave function ^{x, /jL ) 
evolves to (f> w (x, fi), where [i is the characteristic scale of the process. In other words, 
in the leading twist component of the pion wave function the two pion quarks prepare 
themselves for a hard collision by exchanging gluons and increasing their virtualities and 
transverse momenta from the initial scale [i Q ~ Aqcd up to the characteristic scale of the 
process fi < k±, so that this quark pair enters nucleon (nuclei) having already the small 
transverse size r± ~ l/k± (while the smallness of the longitudinal size is ensured by the 
Lorenz contraction). 



3 




Figure 1: Kinematics and notations 
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3 Only the imaginary part of the amplitude, ImT g (from the gluon distribution in the 
target), is shown explicitly in Eqs.(l)-(6), as it is expected to be the main one at high 
energies. For obtaining Im T, the terms it were introduced into all denominators through: 
s — ► s + ie, i.e. £ — > £ — ie [6]. The contribution of the quark distribution to ImT is 
given explicitly in [7]. Let us note also that all integrals over "xi" entering Eqs.(2)-(6) are 
convergent at the ends points X\ — ► or X\ — ► 1. 
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-o = gg (^7^)-^; « = ^ — — j. (7) 

It is seen from the above equations that da/dy is not ~ \(fin(y)\ 2 , but 
depends on the profile of (f> n (x) in a very complicated way. 

4. In this section we present some numerical estimates of the 

cross section, based on the above expressions (l)-(6). Our main purpose was 
to trace the distribution of jets in longitudinal momentum fractions yi,y 2 
depending on the profile of the pion wave function (f) n (x). 

a) For the skewed gluon distribution G^(u,t,fi) of the nucleon at t ~ 
— q\ ~ we used the simple form (as we need it at \u\ > £ only and because 
Gz(u) -> G (u) at \u\ > £): 

G 6 (u,t = 0, li ~ k ± ~2GeV)\ u >z ~ «- a3 (l - w) 5 . (8) 

This form agrees numerically reasonably well with the ordinary, G Q (u,fi ~ 
2Gel / ), and skewed, G^(u,t = 0,/x ~ 2Gel / ), gluon distributions of the 
nucleon calculated in [16] and [17] respectively (in the typical region of the 
E791 experiment: \u\ > £ ~ 1(T 2 ). 

c) As for the pion leading twist wave function, (f> n (x, li), we compare two 
model forms: the asymptotic form, <p® sy (x, fx) = 6x1X2, and the CZ-model 
[18]. The latter has the form: <fi^ z (x, ii Q ~ 0.5 GeV) = 30 #1X2 (^i — x 2 ) 2 , at 
the low normalization point. Being evolved to the characteristic scale of this 
process, li ~ k± ~ 2 GeV, it looks as: <p^ z {x, ll ~ 2 GeV) = 15xix 2 (xi — 

x 2 f + 0.2], see Fig.2. 

The results of these numerical calculations are then compared with the 
E791-data, see Fig.3. 
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X 

Figure 2: Profiles of the pion wave functions: a) ^ z (x,fi ~ 0.5 GeV) = 
30xix 2 (a:i - x 2 ) 2 - dotted line; b) (f>^ z (x,fi ~ 2Gey) = 15xix 2 [0.2 + (xi - x 2 ) 2 } 
- solid line; c) (f)® sy (x) = 6x\x 2 - dashed line. 
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Figure 3: The y-distribution of jets calculated for k± = 2 GeV, E w = 500 GeV and 
with the pion wave functions: 4>^ z (x, fi ~ 2GeV) - solid line, 4>® sy (x) - dashed 
line. The overall normalization is arbitrary, but the relative normalization of two 
curves is as calculated. The data points are from the E791 experiment [1]. 
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It is seen that while the two model pion wave functions are quite different, 
the resulting distributions of jets in longitudinal momenta are similar and, it 
seems, the present experimental accuracy is insufficient to distinguish clearly 
between them. Moreover, even the ratio of the differential cross sections is 
not much different from unity: do asy j da cz ~ 1.2 at y\ = 0.5, and the same 
ratio is ~ 0.7 at y\ = 0.25. Q 

5. The purpose of this section is to compare (in a short form) the 

above results and those obtained by other authors, see [4], [5] and [19]- [20]. 

a) In comparison with Eqs.(l)-(6), the approximations used in [5] cor- 
respond to neglecting the terms S 2 , S3 and £4, and supposing the region of 
small |xi — 3/1 1 dominates the only remaining contribution: / dx\^{x\)lli ~ 
07r(z/i) / dx\Y,\. In this approximation, indeed, ImT g ~ 4> n (y) and da/dy ~ 

\Mv)\ 2 - 

Unfortunately, at the real conditions of the E791 experiment this is a 
poor approximation, and on account of neglected terms in Eqs.(l)-(6) the 
form of da/dy changes qualitatively in comparison with |0^-(|/)| 2 , see Figs. 2 
and 3. f\ 

b) The authors of papers in [4] argued that from the whole set of 31 
Feynman diagrams only a small especially chosen subset dominates ImT g , 
while contributions of all other diagrams will be suppressed by Sudakov form 
factors. From this subset of diagrams they obtained also: ImT g (y) ~ 4>n(y)- 

First, I am unable to understand this last result, as the direct calcu- 
lation of this subset of diagrams according to the rules described above in 
sect. 2 shows that, here also, resulting ImT g (y) depends on (j) n (x) through a 
complicated integrated form, like those in Eqs.(4)-(6) above. 

As for the Sudakov suppression of all other diagrams, I also disagree. The 
line of reasoning in [4] was the following. Consider, for instance, the diagrams 
like those in Fig.2 in [6]. Exchanging the hard gluon, the two pion quarks 
undergo abrupt change of their direction of motion and will tend to radiate 

4 On account of the quark contributions to ImT the resulting curves for da/dy re- 
main very similar to those in Fig. 3, see [7]. For the numerical calculation of ReT see 
[20]. The estimate of ReT can be obtained from the well known formula: ReT(v) ~ 
(■n/2)d\n(v- 1 ImT{v))/d\nv. Because v' 1 ImTiy) ~ zA 3 here, ReT ~ 0.45 ImT and 
|T| 2 ~ 1.2 (ImT) 2 . 

5 The forms of da/dy for (j)^. sv (x) and 4>% z (x, y, ~ 2 GeV) and with the pion energy ten 
times larger, E n = 5TeV, k± = 2 GeV, are shown in Fig. 6 in [6]. It is seen from therein 
that even this does not help much (as the form of the distribution in y of da/dy weakly 
depends on the pion energy), and so the approximation used in [5] remains poor even at 
such energies. 
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gluons in their previous direction of motion. Because the events with such 
radiation are excluded from the data set, what remains will be suppressed 
by the Sudakov form factors. 

Applying this line of reasoning to nearly any hard exclusive amplitude, 
one will conclude that all such amplitudes will be suppressed, as a whole, 
by Sudakov effects. For instance, let us consider the large angle Compton 
amplitude. In the c.m.s. and at high energy, the three proton quarks change 
abruptly their direction of motion, deviating by large angle 9. So, following 
the line of reasoning in [4] , the whole large angle Compton amplitude will be 
suppressed by the Sudakov form factors. 

Really, this is not the case. The reason is that the radiations of three 
collinear proton quarks cancel each other due to colour neutrality of the 
proton. 

c) Unlike [4] and [5], the results from [19]- [20] are very similar to those 
from [6] -[7]. In particular, the analytic expressions for (T g + T q ) in terms of 
corresponding integrals over <p n (x) and the gluon (quark) distributions are 
the same, the only difference is that the signs of ie in denominators of a few 
terms are opposite. For T g for instance (and similarly for T q ), this difference 
can be represented as: 

T g = 5T g (u + ie, M 2 ) + f dx 1 <f> v {x 1 )\f 1 {x i y)g + (v, M 2 ) + f 2 (x, y)g ± (u, M% 

J 



g ± (u, M 2 ) = [' du ^- , (9) 

where 5T g is the same in [6] and [19]-[20], and fi(x, y) and f2{x, y) are definite 
simple functions. The upper sign in Eq.(9) is from [6] (where as pointed out 
above, the terms ie were introduced into all denominators through: v — > 
v + ie), while the lower sign is from [19]- [20]. 

The result for ImT g was obtained in [19]- [20] by direct calculation of 
one-loop Feynman diagrams, and so it can be checked only by a similar 

6 This cancelation becomes ineffective only at the ends points of the quark phase space 
where, for instance, one quark carries nearly all the proton momentum, while other ones are 
wee. The contributions from these regions to the total amplitude are, indeed, suppressed 
by the Sudakov effects. But in any case, the contributions from these end point regions 
give only power corrections to the total amplitude, even ignoring the additional Sudakov 
suppression. 



8 



independent calculation. But there at least one simple argument following 
from general analyticity properties, in favour of the approach used in [6]. 

Let us first consider our amplitude T g (u,M 2 ) in the Euclidean region 
M 2 < 0. There is no discontinuity in M 2 in this region, and the only dis- 
continuity is due to the s-cut. So, ie terms enter unambiguously into all 
denominators as: v = (s — u) — > v + ie, and the terms f\ and f 2 in Eq. (9) 
will be multiplied by the same analytic function g(u, M 2 ): 

[fi(x,y) + f 2 (x,y)]g(u + te,M 2 ) . (10) 

Let us continue now this expression into the Minkowski region M 2 > 0. 
How it may be that after the analytic continuation the functions j\(x,y) 
and f 2 (x,y) will be multiplied by two different functions, as in [19]-[20] (see 
Eq.(9)): 

[h(x,y)g(is + ie,M 2 )+f 2 (x,y)g*(iy + ie,M 2 )} ? (11) 



In conclusion, I would like to emphasize that in spite of some "internal" 
differences between the results from [6]- [7] and [19]- [20], the final curves for 
da/dy are practically the same, and Fig. 3 gives a good representation of both 
answers. So, in any case, the qualitative conclusions agree: under the condi- 
tions of the E791 experiment, the process considered appeared to be weakly 
sensitive to the profile of the pion wave function (f) n (x) and, unfortunately, 
these data can not discriminate between, say, (f)^ sy (x) and (f)^ z (x). 
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